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Abstract. The nature of the radio-quiet X-ray halo around the plerionic SNR G21.5-0.9 is under debate. On 
the basis of spatial and spectral analysis of a large Chandra and XMM-Newton dataset of this source, we have 
developed a self-consistent scenario which explains all the observational features. We found that the halo is 
composed by diffuse extended emission due to dust scattering of X-rays from the plerion, by a bright limb which 
traces particle acceleration in the fast forward shock of the remnant, and by a bright spot (the "North Spur") 
which may be a knot of ejecta in adiabatic expansion. By applying a model of interaction between the PWN, the 
SNR and supernova environment, we argue that G21.5-0.9 progenitor may be of Type IIP or Ib/Ic, and that the 
remnant may be young (200-1000 yr). 

Key words. ISM: supernova remnants; (ISM:) dust, extinction, X-rays: ISM, X-rays: individuals: G21.5-0.9; 
(Stars:) supernovae: general; Radiation mechanisms: non-thermal 



^ 1. Introduction 



The plerionic supernova remnant (SNR) G21.5- 
0.9 has been e x tensiv el y studied in radio (s e e e.g. 
Becker fc Kundul Il976t iBecker fc Szvmkowiakl Il98li 



Fiirst et all ll98St iKassiml Il992t iBock et all 12001 

and references therein) and in the X-ray band (e.g. 
[ Pavelaar et al.ll986tlAsaoka fc Kovamall990tlSlane et alJ 
l2000t IWarwick et all l200lt ISafi-Harb et alJ 12001ft . In 

spite of several efforts, its pulsar remains undetecte d 
llBiggs fc Lvndll99d lLa Palombara fc Mereghettill2002ft . 
IWohier^t^ilTl(]^97]) include it among the non-Crab like 
class of plerions, because of its low frequency spectral 
break, for which a non standard evolution ary path of 
the pu lsar output must be i nvoke d. However. IBock et alJ 
l|200l[) and iBandiera et alJ <)200l[) pointed out that new 
observations, at 94 and 230 GHz respectively, suggest 
a spectral break above 100 GHz, much higher than 
previously thought. 

A set of new and detailed observations in the X-ray 
band performed with XMM-Newton and Chandra have 
apparently raised new and interesting questions about the 
nature of this object, and in particular on the extended 
and diffuse X-ray halo which seems to surround this pie- 
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rion. ISlane et al.l l)200fj(l seem to have been the first to 
detect X-ray emission extending beyond the boundary of 
the radio plerion. The short Chandra calibration obser- 
vation they used prevented a detailed study of the halo, 
and it was not recognized if the emission wa s thermal or 
non-th ermal. Due to the lack of further data, ISlane et alJ 
l|2000|) suggested that the halo may represent the shell 
formed by the interaction of the main blast wave with 
the surrounding medium. They examined archival VLA 
radio data and concluded that the upper limit to the 1 
GHz surface brightn e ss (I ff) is 4 x 10~ 21 W m~ 2 Hz -1 
sr' 1 . IWarwick et all l)200l|) used the XMM-Newton cali- 
bration observation of G21.5-0.9 we also use, and estab- 
lished the non-thermal nature of the X-ray emission of the 
halo. The lack of line emission in the integrated spectra 
of the halo pointed toward a very small ionization time 
(~ 3 x 10 8 cm " 3 s), if non-equ ilibrium of ionization were 
used in the fit. IWarwick et all l)200l[) also noted that the 
size of the X-ray halo exceeds that of the radio PWN 
(pulsar wind nebula) by a factor of 4, a feature which 
is not expected and indeed not observed in any other 
PWN. They also detected a bright spot (named "North 
Spur" ) and some fila ments in the halo. ISafi-Harb et alJ 
l|200lh confirmed the IWarwick et all l|200lh findings and 
considered the halo as an extension of the plerion neb- 
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ula, but they also pointed out that the observed morphol- 
ogy is puzzling and cannot be explained by diffusion mod- 
els. The X-ray m orphology was also briefly discussed by 
iBock et al , who confirmed the previous upper limit 

on Ei GHz, and who suggested multiple events or injection 
epochs, noting, however, that the new radio data do not 
address the question of why the morphology is so peculiar. 
iBandiera fc Bocchinol l|2004h proposed that the X-ray halo 
is due to photons scattered by ISM dust along the line of 
sight, and showed that this model is broadly compatible 
with the observed absorption and radial surface brightness 
profiles. However, their fit to the profiles showed residuals 
due the presence of the (then undetected) weak features 
in the halo, and an error in the profile normalization led 
to the conclusion that an intrinsic halo was require d to 
fit the data in the 5-8 keV. Finally, iBocchind J2005) has 
reported the detection of a weak thermal component in 
the X-ray spectrum of the "North Spur" as seen by both 
Chandra and XMM-Newton. While the metal abundances 
seem to indicate that the North Spur is an ejecta knot, 
the measured X-ray temperature is much lower than the 
expected ejecta temperature behind a reverse shock. 

It is clear that the nature of the X-ray halo of G21.5- 
0.9 is still poorly known. Most of the studies have been 
focused on the central plerion and there has been no sys- 
tematic study of the halo itself. As a consequence, the 
three possible models for the X-ray halo introduced above, 
namely the non-thermal shell, the extension of the plerion 
and the dust-scattering of X-rays from the center, have 
not been properly investigated. In this paper, using new 
X-ray datasets that have been accumulated by Chandra 
and XMM-Newton, we propose a possible scenario which 
explains the halo as a non-thermal shell superposed on a 
dust-scattering profile. We will show that the proposed in- 
terpretation nicely fits all the observational evidence and 
is in agreement with evolutionary models of young PWN 
and SNR. In particular, in Sect. 2 we present the dataset 
we have used, in Sect. 3 we discuss the morphology of 
the X-ray halo, in Sect. 4 we perform spatially resolved 
spectroscopy of the halo and its features, in Sect. 5 we in- 
troduce the dust scattering model which explains part of 
the halo emission, while in Sect. 6 we compare our findings 
with evolutionary models of young PWN-SNRs. 

2. Observations 

G21.5-0.9 was observed as part of the Calibration and 
Perf ormance Verification phase of the XMM-Newton satel- 
lite I .Iansen et alJEoTill) . In particular, the remnant was 
observed both on- axis and off- axis (for the list of G21.5- 
0.9 observations, see Table [TJ- In this work we have used 



Table 1. Cal/PV XMM-Newton Observations of G21.5- 
0.9 



the array of PN CCDs JStruder et a 
arrays of MOS CCDs llTurner et al. 



l200ll) and the two 
2001 !_ at the focus 



of th ree X-ray mirrors of XMM-Newton I Gondoin et alJ 
Il998j) . The nominal bandwidth, angular resolution and 
energy resolution of the instruments are 0.1-15 keV, 15" 
FWHM and E/AE ~ 10, respectively. The original event 
file was screened to eliminate the contribution of soft pro- 



Obs. 


Pointing 


Ti n /T ou t 


(ksec) 


Date 




Location 


PN 


MOS 




OnAxis 


on G21.5-0.9 a 


27/18 


30/25 


7 Apr 2000 


Off Axis 1 


10' S 


28/26 


29/25 


9 Apr 2000 


OffAxis2 


10' W 


28/23 


29/27 


11 Apr 2000 


OffAxis3 


10' N 


29/06 


29/19 


15 Apr 2000 


OffAxis4 


10' E 


29/06 


29/20 


17 Apr 2000 



a At the coordinates 18 fe 33 m 32.6 <i and -10 d 33 m 57 s (J2000) 

tons, from bot h flares and quiescen t emis sion, using the 
recipe given bv lDe Luca fc Molendil l)2004|) . Since the ob- 
servation numbers 3 and 4 were found to be affected by 
unusually large quiescent emission, they were discarded 
in the analysis of the dimmest sources (middle and outer 
halo, see below). All the analysis of XMM data was per- 
formed with the software SAS v6.0. 

G21.5-0.9 was also observed as part of the calibra - 
tion plan for the Chandra satellite ifWeisskoof et al.lll9 96'l. 
There are ~ 70 observations available up to July 2004, 
among which we selected only the 21 observations for 
which G21.5-0.9 was located onto the S3 chip and at an 
off-axis angle less then 5 arcmin. The observation IDs are 
0159, 1433, 1554, 1717, 1769, 1771, 1839, 2873, 3693, 4353, 
5166,1230, 1553, 1716, 1718, 1770, 1838, 1840, 3474,3700, 
4354, and the total exposure t ime is 196 . 5 ks. These are 
the same observations u sed by |Boc^hinc 
includes the set used bv ISafi-Harb et al 



2005). Our set 



2001) plus the 



more recent observations that have been done (they used 
6 observations for a total of 65 ks). Afterward, the data 
were screened for bad grades and for a clean status col- 
umn. All the filtered datasets were merged together using 
the CIAO merge_all task. 

3. X-ray morphology of the halo 

In Figure we show the images of G21.5-0.9 as seen by 
the Chandra ACIS-S camera in two energy bands (0.2- 
2.0 keV and 2.0-10.0 keV, respectively). The figure shows 
the prominent X-ray halo around the bright plerion (the 
latter is marked by radio emission contours) and a com- 
pact bright feature in it at the location 18 ,l 33 m 32.9 ;i and 
-10 d 32 m 49 s (J2000), the so called "North Spur" (NS 
hereafter), located at 80" from the center (~ 1.9 pc at 
5 kpc distance). A search for counterparts of the North 
Spur in other wavelengths performed both with catalog 
browsing (SIMBAD) and image retrieving (Skyview) has 
revealed no obvious counterparts within a radius of 20". 
The circular symmetry of the halo at PA=210°-315° is re- 
markable (here we adopt PA=0 at N, positive clockwise) 
and a bright limb is present. The limb is outlined by the 
white circle in Fig. \I\ which has a radius of 138" (~ 3.3 
pc at a distance of 5 kpc) and a center located 8.7" (0.21 
pc) in the Southeast direction with respect to the peak 
of the PWN. At PA between -45° and 60° we found the 
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NS and bright diffuse emission apparently associated to it. 
This emission is in the form of weak arc-shaped filaments 
which surround the North Spur and connect it to both the 
plerion and the bright limb. In the remaining part of the 
halo (approximately from PA=60°and PA=180°, no other 
features are detected, apart from the star SS397, no limb 
brightening is evident, and the halo declines more rapidly 
toward the background. 

Surface brightness profiles of the halo of G21.5-0.9 
have been derived using the XMM-Newton observations, 
and are shown in Fig. [21 All the profiles are computed 
using a weighted average of the PN, MOS1 and MOS2 
data, where the weights have been derived in each band 
assuming a non-thermal spectrum with a photon index of 
2 and an absorption of 2 x 10 22 cm -2 , which is a good es- 
timate of the halo spectrum. In order to study the profile 
of the "pure" halo, that is without the discrete feature, 
we have selected only data in the PA range 60°-160°. All 
the profiles have been background subtracted using a lo- 
cal background (collected in a large annulus between 280 
and 320 arcsec) and vignetting corrected. Fig. [5] shows 
that G21.5-0.9 has a sharp change in its X-ray profiles 
around 50"(~ 1.2 pc) from the center, where the plerion 
profile (which is narrower than the radio counterpart, as 
expected) suddenly flattens and continues out to ~ 250". 



4. Spectral analysis 

4.1. The halo 

We have performed a spectral analysis of the halo using 
both Chandra and XMM-Newton data. We have gener- 
ated single PN, MOS and ACIS-S3 spectra. The effective 
areas have been averaged, while we have used the stan- 
dard MOS response matrix for the epoch of the obser- 
vations and the on PN matrix. We have defined annulus 
extraction regions, namely region 1 (0"-12"), 2 (12"-24"), 
3 (24"-36") and 4 (36"-51") for the plerion; 5 (51"-94", 
inner halo), 6 (94"-152", middle halo), and 7 (152"-280", 
outer halo) for the halo. Regions containing the North 
Spur and other filamentary structure have been removed. 
We used a power-law emission model to fit the spectra. 
The background has been chosen in an annulus between 
302 and 330 arcsec, excluding out-of-time events for PN. 
The results of spectral fits are shown in Fig. 

The data in region 1-7 are nicely fitted with a power- 
law model, with x 2 /dof = 365/390, 739/395, 338/361, 
338/361, 262/253, 242/205, 181/190, 258/261, respec- 
tively. The derived spectral slope shows a steepening 
in the plerion regi o ns 1—4 whic h was already m e asure d 
by IWarwick et all l)200lj) and ISa.fj-Ha.rh et, ail feOOlh . 
However, in the halo regions 5-7, no significant spectral 
steepening is observed and all these regions are consis- 
tent with a photon index value of ~ 2.3. It is noteworthy 
that the derived value of the absorbing column density 
is not constant in all the regions, as would be expected. 
The absorption is maximum in the central core and de- 
creases monotonically towards the exterior parts. The dif- 
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Fig. 1. Chandra ACIS-S images of G21. 5-0.9 in the 0.2-2 
keV (top panel) and in the 2-10 keV band (bottom panel). 
The color has been chosen to emphasize the weak halo 
emission. Black contours represent 22.3 GHz contours of 
the large-scale diffuse emission of G21 .5-0.9 at 10. 50 an d 



; arcsec HEBW) from lFurst et al.l 



100 mJy/beam 
The position of the North Spur is indicated. The white 
circle fits the outer limb brightened emission of the X-ray 
halo, it has a radius of 138 arcsec (~ 3.3 pc at 5 kpc) and 
a center at 8.7" from the peak of the PWN. 



ference in iVu between the center and the halo periphery 
is AN-r ~ 10 22 cm~ 2 , and there is moderate evidence for 
flattening of the JVh decreasing trend in the halo regions. 
Such a high difference between the core and the halo can 
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Fig. 2. X-ray surface brightness profile of G21.5-0.9 as 
seen by XMM-Newton. For comparison, we also plot the 
profile for a bright point source (LMC X-3). Moreover, we 
plot t he profile of the source as observed bv iFiirst et alJ 
(1988) at 22.3 GHz. 



hardly be intrinsic, corresponding to a line of sight of 3 pc 
of absorbing material at 1000 atom cm~ 3 . 

Moreover, we have also selected a region to study 
in detail: the bright limb of the halo between 115"and 
138" (PA=199°-318°). Spectral fits in these regions were 
performed to test the presence of an additional thermal 
emission due to ISM (interstellar medium) heated by the 
forward shock, using a combination of power-law and 
the mekal model in XSPEC vll.l l|Mewe et all Il985h 
with standar d abundances. W e have also tried the model 
"SRCUT" of lRevnoldsl l)l99S|) which represents the emis- 
sion of electrons accelerated at the strong shock of a SNR 
shell. The latter model was used to tes t if the X-ray hal o 
may be a non-th ermal shell like SN1006 l(Dver et al.l2004h . 
G34 73-0.5 (e.g. lUchiyama et alJl2003h . RX J0852.0-4622 
(e.g. ISlane et alJl200l|) . and other young SNRs. The ther- 
mal component is not detected in the Chandra ACIS-S3 
spectrum, and we derived an upper limit of 10 10 cm~ 5 to 
the emission measure of 1-10 keV plasma, corresponding 
to an upper limit of 0.65 cm~ 3 for the post-shock density 
and an emitting mass < 0.045 M© 1 . In the XMM-Newton 
EPIC spectrum there is a marginal detection of an excess 
below 1.2 keV which requires a thermal component with 
kT = 0.2 — 0.7. Since this is not confirmed in the Chandra 
spectrum of the same region, and since the thermal com- 
ponent does not fit the spectrum in the 0.5-0.8 keV band, 
we do not consider it real. 

As for the fit to the "SRCUT" model, we do not know 
the radio flux at 1 GHz and slope of the radio spectrum (a, 
the energy index), since the halo has not been detected yet 
in radio (for a discussion on the parameters of the SRCUT 



1 Here and in the following we assum e a distance of 5 kpc 
JPavelaar et al]ll98& fBock et al.ll200lTl . The density and the 
mass scale as d" 1 ^ 2 and d 5 ^ 2 , respectively. We have also as- 
sumed the swept-up mass is located in a thin shell. 
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Fig. 3. Results of spectral fitting of G21.5-0.9 plerion and 
halo data with a power-law model modified by interstellar 
absorption. We show the best-fit values and 90% confi- 
dence level uncertainties for the absorbing column density 
and the photon spectral index. 
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0.010 
Normalization (Jy) 



0.100 



Fig. 4. 68%, 90% and 99% confidence level x 2 contours on 
the spectral break vs. radio flux parameters of the SRCUT 
model obtained in a fit of the X-ray spectrum of the bright 
limb in Fig. ^ 
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model see iDver et al-lboO^ . Therefore, we constrained a 
in the range 0.3-0.6, which is typically observed in other 
non-thermal shells, and we left the radio flux and spectral 
break location free to vary. In this way, we may see if 
the extrapolation of the X-ray spectrum back to the radio 
regime according to the acceleration model is consistent 
with the upper limit of ISlane et a l. (2000). The fits of 
SRCUT model to the bright limb are as good as the power- 
law fits, and indicate that the upper-limit, after rescaling 
for the different source regions (the upper-limit is 0.08 Jy 
in the rim region), is ~ 2 — 3 times above the expected 
radio fluxes derived with the spectral fits (see e.g. Fig.rj}. 
The value obtained for the location of the spectral break 
is 2 — 9 x I0 17 Hz, corresponding to a maximum energy of 
accelerated particles of 30 — 85(B/10/iG) TeV, similar to 
what is found in other non-thermal SNR shells. 

4.2. The North Spur 

Th e spectrum of t he North Spur was already studied 
by iBocchind l|2005|) . who reported the presence of an 
additional thermal component. In this work, we test if 
the additional thermal component is affected by Non- 
Equilibrium Ionization (NEI), which is expected to be 
present in the spectra of ejecta, circumstellar and in- 
terstellar material in the SNR. With this aim, we fit- 
ted the North Spur ACIS and EPIC data simultaneously 
with a combination of the power-law model, the VMEKAL 
model of the X-ray emis sion of an opticall y thin plasma 
in ionization equilibrium l(Mewe et alJIlQs"^ . and the con- 
stant te mperature and single ioni zation time NEI emission 
model of iBorkowski et alJ l|200lj) . When fitting the pow- 
erlaw+vmekal and the powerlaw+VNEI combination, we 
fixed the interstellar absorption value to Ah = 2.15 x 10 22 
cm -2 , which is similar to the best-fit value obtained when 
this parameter is left free to vary. The results are summa- 
rized in Tabic [3 while the ACIS and EPIC spectra along 
with their NEI best-fit model are shown in Fig. El The 
X-ray spectrum (especially the ACIS-S3 spectrum) shows 
signs of the presence of two bright emission lines, namely 
Mg XI at 1.34 keV and Si XIII at 1.86 keV (Fig®. 

Unfortunately, the combination of the power-law and 
NEI models gives two local minima, one which repre- 
sents the equilibrium situation already found with mekal 
fit (x 2 /dof = 626.7/583, hereafter Ml), and one which 
represents a plasma strongly affected by NEI conditions 
(x 2 /dof = 628.8/583, hereafter M2). Both results are pre- 
sented in Table El and displayed in Fig. El 

The minimum Ml, which has a slightly lower \ 2 value 
then M2, yields parameters similar to the ones found by 
iBocchinol l)2005h . On the other hand, the minimum M2 
gives an ionization time between 7 x 10 9 and 2 x 10 10 
cm~ 3 s and a temperature in the range 0.2-0.4 keV (at 
the 90% confidence level), a factor of two higher than Ml 
(Fig. E| top panel). 

The NEI results at minimum M2 lead to an emitting 
plasma density of 5(3 — 8) cm~ 3 and mass of 0.23(0.17 — 



X — ray spectrum of the North Spur 



EPIC PN (top), ACIS — S3 (middle), EPIC MOS (bottom) 
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Energy (keV) 



Fig. 5. Chandra ACIS-S3 and XMM-Newton PN and 
MOS spectra of the North-Spur, along with the best-fit 
NEI model. We also show the contribution of thermal 
and non-thermal component to the ACIS spectrum (dot- 
dashed and dashed lines, respectively). 

0.37) M Q , if the line of sight extension of the emitting 
plasma is equal to the chord intersecting the sphere of the 
SNR shell (~ 6.5 pc) At the minimum Ml the correspond- 
ing values are ~ 10 cm -3 and ~ 0.5 M Q , respectively. It is 
noteworthy that the derived age of the North Spur com- 
patible with the M2 minimum is ~ 100 yr, which is signif- 
icantly lower than the corresponding age at Ml ( > 1000 
yr, Fig. E| bottom panel). 

The metal abundances of Mg and Si, as measured by 
the minimum Ml are 0.6-3 times the solar value for Mg 
and 2-20 for Si, but for M2 the abundances are con- 
sistent with the solar values for both elements (Fig. 0). 
The measured thermal flux corresponds to a luminosity of 
- 4 x 10 34 erg s^ 1 in the 0.5-2.0 keV band. This lumi- 
nosity is a bit high compared to an eje cta knot of Cas A, 
for which from the results published by Laming & Hwang 
l|2003|) we have computed L max ~ 10 34 erg s _1 . However, 
our measurement of the unabsorbed thermal flux is af- 
fected by large uncertainty, which can increase or decrease 
the luminosity by a factor of 10. The total absorbed flux 
of the North Spur in the 0.5-2.0 keV band is 3.2 x 10" 13 
erg cm -2 s _1 , of which ~ 25% is due to the thermal com- 
ponent. 

5. A dust scattering model for the "pure" halo 

In this section, we investigate whether the diffuse emis- 
sion of the G21.5-0.9 X-ray halo may be due to dust 
scattering of X-ray photons from the plerion. There are 
reasons to suspect that dust scattering contributes appre- 
ciably to the G21.5-0.9 X-ray profile. One is the large 
absorption column density toward this source. Dust scat- 
tering X-ray halos are typically found around other heav- 
ily absorbed sources (see e.g. IPredehl fc Schmittl Il 995: 
ISmith fe Dwek1ll998t iNaease et all l200lt IVaughan et alJ 
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Table 2. Results obtained by a ACIS-EPIC joint spectral fit of the North-Spur X-ray emission. We fixed the interstellar 
absorption at iV H = 2.15 x 10 22 cm' 2 . 

Model" 7 Norm kT t thermal flux" X T Jdof 

10 22 cm" 2 photon index ph cm -2 s _1 keV -1 keV cm -3 s erg cm -2 s _1 

~TZ 2.45 ±0.05 2.7 ±0.2 x 10" 4 - - - 876.0/589 

PL±VMEKAL 2.18 ±0.04 2.0 ± 0.2 x 10 -4 0.13 ±0.06 - 2.3 x 10" 11 632.8/584 

PL±VNEI (Ml) 2.15 ±0.04 1.8±0.2xl0 -4 0.17(0.15-0.21) 7 x 10 11 1.5 x 10 -11 626.7/583 

PL±VNEI (M2) 2.15 ±0.04 1.9±0.2xl0~ 4 0.30(0.20 - 0.37) 1 x 10 10 1.7 x 10~ n 628.8/583 



a PL=power-law 

6 Unabsorbed flux in the 0.5-2.0 keV band due to the thermal component only. 




0.2 0.3 0.4 0.5 



Temperature (keV) 

Post-shock density (cm 3 ) 
3 5 10 20 



■ - 


i 


< 1 \ \ "\ 1 : 
\ M1 (1075) )\ r ~~\~' '■ 










5 M2 (1 .079) 1 


i 



' ' ' 'i'o' 12 ^ 

Normalization (cm 5 ) 



Fig. 6. Ranges of kT and ionization time (top panel), 
emission measure and normalization time (bottom panel) 
allowed by the fit to North Spur data with a combina- 
tion of power-law and NEI emission model. In the bottom 
panel, we also show the range of post-shock density and 
the range of ages which are compatible with the EM and 
r values derived from the fit. 

l2004|b IPredehl & Schmittl (|l995|) have shown that there 
is a good correlation between the optical depth for dust 
scattering and the absorption column density. Using their 
empirical formula, for G21.5-0.9 one expects TikcV — 1 
(where t\ keV is the scattering optical depth a 1 keV) . For 
the sake of illustration, n koV = 1 implies a halo frac- 
tional flux of about 50%, 20%, 7% and 3%, respectively 
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Fig. 7. Mg and Si abundances obtained by fitting with 
a powerlaw and VNEI model. The locations of the two 
minima (Ml and M2) are marked by arrows, and the % 2 
is reported. 

for the energy bands 1.3-1.7 keV, 2-3 keV, 3-5 keV and 
5-8 keV (these will be the reference energy bands, used in 
our analysis). However, there is a large scatter about that 
correlation, so that it cannot be taken as a safe method 
to evaluate r\ k c v m individual sources. 

Another effect that could be also ascribed to dust scat- 
tering is the systematic decrease of the absorption column 
density with the distance from the source center (see Fig. 
|3J), as derived from X-ray spectral analysis. If it was a 
true column density variation, it would require an im- 
probable distribution of the foreground matter. Instead, 
a lower measured absorption in the outer regions can be 
justified as an artifact of the spectral analysis, in the pres- 
ence of dust scattering, due to the fact that hard X-rays 
are scattered at lower angles than soft X-rays. 

We present here an analy sis along the lines of that by 
iBandiera fc Bocchino! l)2004|b but carried out with much 
better accuracy. First, we start from data of better quality, 
because they are based on a longer integration time, and 
with more careful selection of the directions along which 
the contamination from the North Spur and the shell is the 
lowest. Then, we have corrected for some small bias in the 
profile normalization, which was affecting the analysis of 
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iBandiera fc Bocchinol l|2004j) . Finally, we have performed 
a much better analysis on the parameter space. 

In short, the procedure is as follows. We assume that 
the dust halo is negligible in the 5-8 keV band (we have 
already seen that, for TikoV = 1, the halo fractional flux 
in this band is only 3%). Therefore, the radial profile in 
the 5-8 keV band is well approximated by the convolution 
of the intrinsic source profile with the instrumental Point 
Spread Function (PSF). Let us assume that the PSF is en- 
ergy independent (this is a good approximation, for EPIC 
onboard XMM-Newton) and also that the source intrinsic 
radial profile does not change much with the energy (this 
assumption will be verified "a posteriori" ). In this case, the 
5-8 keV radial profile should reproduce rather well profiles 
in other bands, except for the different flux normalizations 
and for the further scattering component, which becomes 
more prominent in softer energy bands. For a given model 
of the scattering halo, the scattering component is thus ob- 
tained as its convolution with the PSF. Finally, the sum 
of the intrinsic and scattering components must be com- 
pared with the profile measured in a given energy band 
(we shall use 2-3 keV). A least x 2 analysis allows one 
to constrain the parameters of the halo model. This is a 
rather time-consuming procedure, because for each set of 
parameters we are required to compute one convolution of 
2-dimensional maps. 

It is well known that, as long as the Rayleigh-Gans 
approximation is valid, the scattering optical depth scales 
as E~ 2 (where E is the photon energy), while the angu lar 
size of the halo scales with E~ x (see e.g. |Predehllll998l for 
a discussion of these scaling laws). For the angular size, 
a formula similar to that of classical diffraction applies, 
namely scattering angles are of the order of the ratio be- 
tween photon wavelength and dust grain size. Using this 
scaling law, we can then simulate profiles for the other 
two energy ranges, 1.3-1.7 keV and 3-5 keV, and com- 
pare them with observations. 

However, it should be clear that even in this sophis- 
ticated analysis we can only check for self-consistency of 
a scattering halo model, while we cannot prove for sure 
the dust-scattering nature of an X-ray halo. In princi- 
ple, both spatial and spectral properties of a source may 
conspire to mimic shape and scaling of a dust-scattering 
halo. However, the required conditions would be highly 
unlikely, and this is the reason why we finally conclude 
that the dust-scattering halo hypothesis is confirmed to a 
high confidence level. 

The relevant parameters to model a dust-scattering 
halo are the following: the optical depth for scattering 
(at the reference energy, say rikcv); the dust distribution 
along the line of sight (here we shall use a uniform dis- 
tribution); the size distribution of grains, usually approx- 
imated by a power-law with an upper cutoff (oc a~ q for 
a < a max ; typical values are a max = 0.17 /im and q = 3.9, 
Predehl and Schmitt 1995). 

For the scattering halo, we use a simplified, analytic 
model, that will be described in a forthcoming paper, to- 
gether with details of the procedure of halo modeling and 



CO 

w 

CD 

c 10" 

-t-> 

jl 

M 

£ 10" 

CD 
O 
CO 

t 10" 
=s 

tn 

-a 

v 10" 



CO 

u 10" 
o 




50 100 150 200 250 
Distance from center (arcsec) 



0.4 



% 0.3 
CO 

I 

CV! 



■0.2 



0.1 



iii! 


i i i | , i , | , i 

uniform dust 


distr. 




q=3.9 






errors meas 


values 




y- "7416^0. 143)--^' 






\ V \ x ^^^^ 




- , , , 1 


, , , l , , , l , 


i , , , - 



200 



400 600 

0=„»i(arcsec) 

scal^ ' 



800 



Fig. 8. Top panel: G21.5-0.9 X-ray profiles with best-fit 
dust scattering model (crosses 1.3-1.7 keV, asterisks 2-3 
keV, diamonds 3-5 keV, triangles 5-8 keV). The derived 
dust parameters and their uncertainties are shown in the 
bottom panel. 

subtraction in individual sources. We outline here only 
some basic features of this model and of halo models in 
general (under Rayleigh-Gans, and single-scattering ap- 
proximations). 

If the space distribution of dust is uniform along the 
line of sight 2 then the halo profile is oc 8" 1 , for 9 <C sca \, 
where 



0scal(£) 



3200 arcsec 



(£/lkeV)(a max /0.17^m) : 



(1) 



2 We have tried to relax this hypothesis by considering a 
thin slab of scattering dust at a given distance. The chi-square 
analysis shows that the "uniform" case has to be preferred. 
More elaborate and complex scenarios are not justified by the 
present data quality. 
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Fig. 9. Image of G21. 5-0.9 in the 2-8 keV with dust scat- 
tering contribution subtracted. 



The range of radial distances over which we see the 
halo in G21.5-0.9 is about 100-300 arcsec. Therefore, from 
Eq. n it follows that we are in the 9 <C 9 sca i regime. A 
consequence of this (as we have found by analyzing nu- 
merically the parameter space) is that the best fit halo 
does not constrain at all the value of q: then, we have as- 
sumed for this parameter the average value 3.9, as given 
by Predehl and Schmitt (1995). 

For the minimum \ 2 fit, we have used errors propor- 
tional to the measured values, because the radial profile 
within ~ 100 arcsec is much more uncertain than one 
would infer from the (very small) statistical error (here 
the main uncertainty derives from the fact that the in- 
trinsic source profile slightly changes with energy). The 
best-fit we got is shown in Fig. |S| along with the de- 
rived t (in the 2-3 keV band) and 9 sca \. We thus esti- 
mate TikcV = 5.68r2_3keV — 0.80, with a la uncertainty 
of about 20%. 

In Fig. Elwe show an image of G21.5-0.9 in the band 
2-8 keV with the best-fit dust scattering model profile sub- 
tracted. The emission which remains after the halo sub- 
traction is the residual flux from the plerion (which is not 
expected to be properly modeled by a circularly symmet- 
ric model), the North Spur, its filaments and the shell. 
In particular, the shell is clearly visible from PA=180° to 
300°, while it is not visible at PA=90°-180°. 

In a separate paper, we will report more details on the 
dust scattering model and the fit to the G21.5-0.9 X-ray 
halo, including a thin slab dust distribution and a discus- 
sion of the dust properties we found, in the framework of 
general dust properties as derived from studies of X-ray 
halos in other sources. 



6. A self-consistent model for the X-ray halo and 
its features 

6.1. Establishing the evolutionary stage of G21.5-0.9 

As we have seen, the bright limb in Fig. rj] and El suggests 
that we have detected the forward shock of G21.5-0.9 ex- 
panding into the environment of the SNR. If this is the 
case, we may use evolutionary models of young SNR to 
infer the age of G21.5-0.9. First, we argue that the rem- 
nant is not yet in the Sedov- Taylor (ST) phase, because 
the swept-up mass by the shock front in the bright limb, as 
derived in Sect. 4, is very low (< 0.045 Mg). Even taking 
into account the mass in the North Spur and the other fil- 
aments not analyzed, the total swept-up mass in the halo 
is less than 0.5 M . 

iGhevalierl ll2005h has discussed the interaction of a 
young core-collapse SNR with its environment, taking into 
account different supernova types, namely SN1987A-like 
class, IIP, IIL/b, and Ib/Ic. Because of the low mass ob- 
served in the halo, we may safely discard a dense environ- 
ment for G21.5-0.9 like the RSG wind usually found for 
Type IIL/b. In fact, it can be shown that for an observed 
radius of 3.3 pc (Sect. 3), the remnant of a Type IIL/b 
with a dense RSG wind should have swept-up most of the 
total ejecta mass and a few solar masses of circumstellar 
medium. 

Type IIP and Ib/Ic ar e interesting possibilities. For 
Type IIP, IGhevalierl §005) notes that the low mass loss 
during the RSG phase would result in a small region (r < 1 
pc) of dense wind surrounded by a more diffuse extended 
bubble created during the main sequence phase. The in- 
teraction with the resulting double layer wind structure 
has not been modeled in detail but swept-up masses are 
of the order of 0.1 M , so in agreement with observa- 
tions. For Type Ib/Ic, it is expected that they undergo a 
Wolf-Rayet star phase with high mass-loss and fast wind 
which sweeps the earlier RSG bubble in a shell at sev- 
eral parsecs from the center. The evolution of the SNR 
in this medium ma y also be complica ted and requires nu- 
merical simulation l|Dwarkadasll200ll) . However, if we as- 
sume that the star has been a WR object long enough 
to produce an extended wind co mponent, we may use the 
interaction model developed bv IChevalierl {2005) for the 
circumstellar interaction of an RSG wind, but with the 
parameters appropriate for a WR wind (-B51 = 1, M c j = 4 
M , M = 3 x 10 5 M yr" 1 , wind velocity -y w = 1000 km 
s" 1 , D = M/4tto w = 1.5 x 10 12 g cm" 1 ). The model as- 
sumes that the progenitor star had a radiative envelope 
and treats the interaction in the thin shell approximation. 
Although the progenitor star of G21.5-0.9 may not have 
had a radiative envelope, the steep outer power law with 
a relatively flat central density distribution is probably a 
reasonable approximation to the density profile. We as- 
sume that the interaction shell is still in the outer steep 
power law part of the supernova density profile; this as- 
sumption can be verified for the parameters we find for 
G21.5-0.9. In this case, we derive a CSM swept-up mass 
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Fig. 10. The CSM interaction model of IChevalieil l|2005|) 

adapted to the G21.5-0.9 case. The two segments are the 
loci of allowed wind parameter and SNR age values calcu- 
lated for 4 and 10 M ejecta masses and an upper-limit of 
0.1 Mq for the CSM swept-up in a shell, in agreement with 
the results obtained in Sect. 4.2 for the bright halo. The 
wind parameters for three kinds of stars are also shown, a 
RSG wind (which turns out to be incompatible with the 
data for any reasonable value of the system parameters) , 
a WR star wind (one of the viable interpretations for the 
G21.5-0.9 progenitor), and a 04V normal star, shown only 
as a reference for the minimum value the wind parameter 
may have. 

of 0.1 Mq, an ejecta swept-up mass of 0.2 M©, which 
are again roughly in agreement with observed mass in the 
North Spur and the upper-limit of the shell. The derived 
shell velocity is ~ 7000 km s -1 , and the age is ~ 290 yr. 
For any reasonable set of wind parameters, the upper- limit 
on the swept-up mass derived in Sect. 4.1 implies an age 
lower then ~ 450 yr. This is shown in Fig. ^| where al- 
lowed SNR ages are plotted for different wind parameters 
and ejecta masses according to the self-similar s olution of 
the CSM interaction model of IChevalierl l|2005h . adapted 
to the observed size of G21.5-0.9, and considering only a 
swept-up CSM mass below 0.1 Mq, in agreement with the 
observations. 

Therefore, the spectral results for the shell and the 
North Spur, when combined with reasonable scenarios of 
the SN environment, seem to indicate in any case an 
early evolutionary stage for the SNR and a progenitor 
without extended dense winds. We may now see if this 
is also compatible with the size of PWN and some ap- 
propri ate model for PWN-SNR interaction. Several au- 
thors llRevnolds & Chevaherlll984l IChevalier & Franssor] 
1992: van der Swaluw et alJl200l|) give analytical expres- 
sions which relate the radius of the PWN shock and the 
SNR blastwave. We have applied these equations to the 
case of G21.5-0.9 using a SNR shock radius of 3.3 parsec 
and a PWN shock radius of 1.2 parsec, under the assump- 
tion of a uniform ejecta density and constant spin-down 
luminosity. We typically get E pv!n = L t ~ 10 49 erg and 



Table 3. A comparison of G21.5-0.9 linear size and 
equipartition energy with a set of similar plerion with 
known age. From the comparison, one would guess a 
G21.5-0.9 age in the range 800-1600 yr. 



SNR 0540-69 Kes75 Gil. 2 G54.1 G21. 5-0.9 

R (pc) 0~9 O 09 L2 L2 

Age 800 1000 1600 1500 ? 

E min 5 x 10 47 10 48 3 x 10 46 8 x 10 46 4 x 10 47 



an age of ~ 500 yr from these models. While the age is 
in agreement with the CSM interaction model, the cor- 
responding (constant) spin-down luminosity is very high 
especially considering that the X-ray luminosity is consid- 
erably less than that observed from the brightest PWNe, 
so this seems to suggest that significant spin-down has 
oc curred at the earl y stage of the PWN evolution. 

IChevalierl l|2005|) has developed a model of PWN in- 
teraction with SN ejecta, considering power-law density 
profile for the ejecta and the pulsar luminosity decay. In 
the framework of his model, it is convenient to compare 
the PWN internal energy with the pulsar luminosity and 
the shell kinetic energy. The equipartition energy of the 
nebula can be found in the following way 

4 9 

E min = 3.2 x 10 47 ( J*-) * ( **™) 7 (—%-) ^ (2) 
V3-6Jy/ V !- 2 P c / V540GHz/ w 

where Lb is the spectral luminosity a t the break fre- 
quen cy fb (values for G21.5-0.9 are from iBandiera et all 
200lh . the particle energy indexes before and after the 
break are r»i = 1 and p-2 — 3, respectively (adapted from 
IChevalierl 12001 . The nebula internal energy is usually 
within a factor of a few of E min , say E int — 10 48 erg. 
ISafi-Harb etafl fcOOlh found £ = 3 x 10 37 erg s" 1 , based 
on the E — Lx relation, implying Et = 4.7 x 10 47 isoo erg, 
where isoo is the age in units of 500 yr. The kinetic energy 
in the PWN is ~ 10 49 erg, which is determined by the 
supernova model and the age. This set of energies is con- 
sistent with the model if the ratio of the age and the initial 
spin-down time scale tjr is in the range 1-10 (see Fig. 1 in 
IChevalierl2005V thus suggesting that significant spindown 
can occur at the early stage of the PWN evolution. 

Finally, we should note that the comparison of the 
equipartition energy and linear dimension of G21.5-0.9 
with other PWNe, summarized in Table [3] seems to sug- 
gest an age of 800-1600 yr, older then what we found 
on the basis of the circumstellar interaction model. One 
way of having the PWN expand more rapidly is to have 
a lower density supernova, but this generally requires a 
lower mass or higher energy, which also reduce the age 
from the circumstellar interaction model. 

6.2. On the nature of "North Spur" 

In this section, we review a physical interpretation for the 
North Spur which is compatible with the observational 
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results. Fig^shows that the spur is located at 35" (0.8 pc) 
outside the plerion rim, with a low surface brightness in 
between. This is at odds with an interpretation in terms of 
material swept-up by the plerion, because in that case we 
expect the thermal emission adjacent or i nside the plerion 
non-t hermal emission (as in, e.g., 3C58, iBocchino et alJ 

Therefore, a possible valid interpretation may be that 
the North Spur is an ejecta clump hit by the reverse shock. 
In this case, the position of the forward shock (Rf s ), as 
traced by the bright limb in the halo, and the contact dis- 
continuity (i? c )j as traced by the emission from the North 
Spur, can be compared with self-similar solutions of inter- 
action of e j ecta w ith an external medium as worked out by 
IChevalierl l|l982() . For a model in which ejecta with steep 
power-law outer density distribution (n > 5, where n is 
the power-law index) are interacting with the free stel- 
lar wind of the massive progenitor (s = 2, where s is the 
power-law index of stationary medium, usually mass-loss 
from progenitor star), the forward shock is expected to 
be at < 1.3 — 1.4i? c , unless it has already entered the 
flat part of ejecta density distribution. Unfortunately, be- 
cause of projection effects, the exact location of the North 
Spur inside the remnant is not exactly known, but it is 
between 1.0i? c and 1.7i? c , so in general agreement with 
expectation. As already noted in Sect. 4.2, the density, 
mass and luminosity estimates for the North Spur are not 
unusual for ejecta in young SNRs. However, the measured 
temperature is lower then the one expected for emitting 
X-ray ejecta and usually observed in oth er young SNRs 
originating from a Type II SN (e.g. Cas A. lGotthelf eTaD 
200l|) . thus casting some doubt on this interpretation. 

A different explanation for the North Spur which takes 
into account the low X-ray temperature may be the fol- 
lowing. We have seen the a progenitor SN of Type IIP 
is in agre e ment with the observationally derived masses. 
IChevalierl l)2005|) show that in case of a Type IIP SN, the 
circumstellar interaction should occur early (radius < 1 
pc), followed by adiabatic expansion of the reverse shocked 
ejecta. This adiabatic expansion could give rise to the rela- 
tively cool emission which is observed. The mass is not ex- 
pected to be high and it is in agreement with the observed 
value of ~ 0.2 M & . Other implications of this scenario are 
a young age and near solar abundances (the envelope of 
the the Type IIP SN), which are both in agreement with 
the age and abundances derived at the minimum M2 (Fig. 
EJ), and the dynamical age of the shell derived in the pre- 
vious section. 

We also note that the data indicate that the North 
Spur has an intrinsic non-thermal tail in the spectrum 
that is modeled with a power-law (7 = 2.15, Table El- One 
explanation for the non-thermal emission is synchrotron 
radiation from particles accelerated in a shock. 

7. Conclusions 

We presented an extensive analysis of Chandra and XMM- 
Newton X-ray data of the radio-quiet halo around the 



plerion G21.5-0.9. We included in our analysis all the 
public observations of the source available up to now. 
We showed that the halo morphology is composed of two 
components: diffuse emission and some bright knots and 
filaments. We confirmed the detection of X-ray thermal 
emission in the brightest kn ot, the "North Spur", already 
detected bv lBocchind (|2005fl . and we presented a more de- 
tailed spectral analysis of this object which included Non- 
Equilibrium Ionization, and which yielded a very young 
age (100-1000 yr) and abundances compatible with solar 
values. We have presented various interpretations for the 
origin of this knot, and we argue that it can be due to 
ejecta interaction with the H envelope of a Type IIP SN. 

We have also detected a bright limb in the east part 
of the halo, located at 3.3 pc from the center. The limb is 
dominated by non-thermal X-ray emission, probably due 
to particle acceleration at the fast forward shock, while 
the upper limits for the post-shock density and emitting 
mass are 0.65 cm -3 and 0.045 M Q in this region. 

We showed that the diffuse emission from the halo is 
due to dust scattering of X-rays from the plerion, and we 
explored which range of dust parameters is compatible 
with observations. 

We argued that our data are inconsistent with an ex- 
planation of the halo in terms of a plerion extension, as 
suggested by previous works. We exclude the possibility 
that the system is in Sedov- Taylor stage on the grounds 
that the swept-up mass is very low. Moreover, by apply- 
ing a model of CSM interaction of a young SNR shock 
wave to the bright limb data, we argue that G21.5-0.9 is 
in a very early evolutionary stage, with an age between 
250 and 500 yr. A comparison of the plerion linear size 
and equipartition energy with the corresponding values of 
plerions with known age suggest a longer age (800-1600), 
while the comparison with a model of SNR-PWN interac- 
tion for young remnant allows age as low as 500 yr and 
indicates that spindown has already occurred. Putting to- 
gether all the estimates, it seems that a very reasonable 
range for the age of G21.5-0.9 is 200-1000 yr. The lack of 
an historical supernova associated with this object is not 
a surprise given the high extinction in this direction. 

If our conclusions are correct, the PWN may be ex- 
panding at a few times 1000 km s _1 , and a comparison 
of X-ray archive images at a baseline of 5 yr may lead to 
the direct detection of PWN expansion. Moreover, addi- 
tional X-ray observations would be required to study in 
more detail the thermal emission both from the forward 
shock (if any) and in the North Spur. Deeper radio ob- 
servations of the North Spur and the bright limb which 
lead to detection of these objects would shed light on the 
nature of their non-thermal emission. The required radio 
sensitivity at 1 GHz to detect the limb is 2-3 times below 
the current upper- limit. As for the halo X-ray diffuse ex- 
tended component, if it is indeed due to dust scattering, 
it should not have any radio counterpart. The detection 
of the pulsar would be of great value to further constrain 
the evolutionary stage of the system. 
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Note added after acceptance: In April 2005, dur- 
ing the referee review of thi s manuscript, a paper by 
Matheso n fe San- Harb (2005) appeared, reporting the 
results of a spectral analysis based on a Chandra 
da taset larger then the one u sed by us. The results 
of iMatheson fc Safi-Harhl l(2005h are consistent with our 
ones, in particular the trend in the interstellar absorption 
and power-law index, and the presence of limb brightening 
in the east quadrant. 
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